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تاریخ دریافت: ۱۳۹۵/2/۲۶ تاریخ تصویب: ۱۳۹۵/۱۲/۲۵ 


| ستفاده از مدل‌های هیدرولوژیکی یکی از روش‌های رایج در تجزیه‌وتحلیل کمی حوضه‌های آبخیز است. 
پیشرفت‌های اخیر در زمینه‌سازی‌های بارش "رواناب که انعطاف بیشتری در حل مسائل و پدیده‌های 
هیدرولوژیکی دارند. آن‌ها را جایگزین مناسبی به‌جای روابط تجربی کرده اسست. به‌طور مرسوم در 
پژوهش‌های صورت گرفته پیشین مبتنی بر بهینه سازی مدل در همان مقیاس محلی موردبررسی قرارگرفته 
که فاقد عمومیت مکانی اتستکت : هدف این مطالعه ارائه پارامترهایی با اندازه معین ات تفش تا در سطح 
گراف سیلاب که در هیدروگراف سیلاب نقش دارند. در مقیاس زیرحوضه‌ای و بزرگ‌مقیاس مورد مقایسه 
راز گرفته است, ور این مطالغه از سل 0۳3۸5 یرای مل ساری رواناب ویر عضو طه‌ها و میتان ووانانب 
خروجی از حوضه استفاده‌شده است؛ برای ساختار مدل از عوامل شماره منحنی 5 برای روش تلفات و 
هیدر و گراف واحد حهت روش انتقال استفاده شد. به‌منظور واسنجی ابتدا بهینه‌سازی پارامترها در 
هریک از زیرحوضه‌ها به صورت جداگانه انجام گرفت و سپس مقدار عمومی پارامترها با به د.ست آوردن 
اندازه ثابتی از پارامترهای حساس که در همه زیر حوضه‌ها پاسخ قابل قبولی ارائه دهنده انجام شسد. 
استفاده از پارامترهای بزرگ‌مقیاس است. ولی ایستایی پارامترهای عمومی بهتر اسشنستت. به‌منظور ارزیابی 
کارایی مدل از شاخص 0101186:-90هآ"استفاده گردید که این شاخص برای کالیبراسیون محلی ۰/۸۵ و 
کالیبراسیون بزرگ‌مقیاس ۰/۱۵ به دست آمد که در دامنه مطلوبی برای شبیه‌سازی قرار دارد. 

کلیدواژه‌ها: مقیاس محلی» بزرگ‌مقیاس» مدل هیدرولوژیکی, 1376-115 505. 
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۱. مقدمه 


یکی از روش‌های رایج کمی در تجزیه‌وتحلیل حوضههای آبخیز بهره‌گیری از مدل‌های هیدرولوژیکی است. 
دسته‌ای از اين مدل‌ها برای شبیه‌سازی فرآیند بارش- رواناب مورداستفاده قرار می‌گیرند. یک مدل هیدرولوژیکی 
بارش -رواناب با ادغام داده‌های سری زمانی؛ پارامترهای مثر. متغیرهای منطقه و قوانین فیزیکی در داخل سیستم 
حو ضه آبخین شبیه سازی میزان رواناب و فرآیندهای دیگر را امکان‌پذیر می سازد. هیدروگراف سیلاب نما شی 
گرافیکی و مهم در تحلیل‌های هیدرولوژیکی می‌باشد. پاسخ‌دهی حوضه در شکل هیدروگراف با مدیریت ارتباط 
مستقیم دارد. در واقى مدیریت حو ضه آبخیز مسر نخواهد بود مگر اینکه خحصو صیات هیدرولوژیک حو ضه برای 
درک سیستم پیش‌بینی شود (بهادرا و همکاران ۲۰۰۸: ۱۰۱۸). 

مدل‌های هیدرولوژیکی متنوعی برای نیازهای مختلف ارائه‌شده‌اند. مدل‌هایی از قبیل مدل آبخیزداری دانشگاه 
کلمبیا ( 611307۷۷( وکاس و وازیلیادس:۲۰۱۶) و مدل ارزیابی آب‌وخاک (8۷۷۸۲۳(سرویس پژوهش کشاورزی ۲ 
ابالاتمتحده ۱۹۹۰ هیدروگراف واحد لحظه‌ای ژتومورفولوژی (01113) ب‌صورت جهانی در حوضه‌های فاقد 
آمار | ستفاده می شود (ابراهیم باتیش و احمد۲۰۱: ۱۳). اين قبیل مدل‌ها برای شبیه سازی فرآیند بارش -رواناب 
بیش از ۶ دهه اخیر مورداستفاده قرار می‌گیرند؛ اما ارائه رواناب برای حوضه‌های فاقد آمار هنوز برای هیدرولوژیست‌ها 
به‌عنوان یک چالش باقی‌مانده است (واسوا"۱۳ ۰ ۲۰). همین امر منجر به تنوع در مدل‌ها شده است. 

مدل3۷]5]-)1۳ یکی از مدل‌های شبیه‌سازی فرآیند بارش - رواناب است. مدل 11۳0-115توسط مرکز 
مهند سی هیدرولوژی ارتش آمریکا ارائه شده ا ست که از قابلیت‌های مهم آن امکان ارتباط با دیگر نرم‌افزارها ازجمله 
سیستم اطلاعات جغرافیایی است. الحاقیه 11015 1۳0-050 ابسیاری از داده‌های مورد نیاز را که نیاز به محاسبه 
دستی دارند. اتوماتیک محاسبه می‌کند. مدل به صورت کاربر پسندی. حوضه را با آیکون‌های مسیرها و اتصال‌ها و... 
موجود در محیط ۳1۳-۳115 برای ورود به مدل آماده می‌کند. این مدل هم قابلیت اجرا به‌صورت یکپارچه و هم 
به‌صورت توزیعی دارد. کابرال موریرا (۲۰۱۳: ۱7۶) به‌منظور مدل‌سازی رواناب در حوضه رودخانه گرانجریو از مدل 
نیمه توزیعی 11۳80-11۷15 برای به د ست آوردن پارامترهای موردنیاز استفاده کرد و مدل حاصل را مدل19110نامید. 
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همچنین آگراوال (۲۰۰۵: ۱۷۰) به‌منظور ارائه مدلی برای ورود داده‌ها به 1۳0-112۷16از مدل توزیعی استفاده کرد و 
مدل حاصل ر ۳0004 نامید. 

از نکات این مدل امکان تغییر ساختار سیستم مدل‌سازی توسط کاربر است. در مدل 1۳0-۳115 ]برای محاسبه 
روش تلفات. انتقال و روندیابی جریان روابط مختلفی وجود دارد که هر رابطه پارامترهای خاص خود را دارد که با 
| ستفاده از روابط تجربی محا سبه می شود. این مدل می‌تواند در مناطق فاقد آمار و در مناطقی که نظارت روی داده‌ها 
محدود است میزان رواناب را برآورد کند (ابراهيم باتیش و احمده۲۰۱۹: ۱۱). 

دقت برآورد رواناب با مدل‌های هیدرولوژیکی به‌دقت برآورد پارامترهای مدل و متغیرهای مکانی و زمانی بستگی 
دارد (کاراپوا" و همکاران۲۰۱۲: ٩۵٩‏ لو" و همکاران۲۰۱۳:۱۱۰؛ بویج*۲۰۰۳: .)۲۵٩۱‏ درواقع می‌توان گفت میزان 
موفقیت در تعیین دقت شبیه‌سازی وابسته به‌دقت پارامترهای ورودی و شناخت کاربر از ساختار مدل است 
(زانگ ۲۰۱۳: 1۵؛ بورمن »۲۰۰۳ ۱۷۶؛ باتس ۲۰۰۶: ۲۷۰). راه‌حل بهینه این است که کالیبره کردن داده‌های بارش 
با توجه به تغییرات زمانی و مکانی باشد. به دلیل نقش کلیدی بارش اندازه‌گیری‌های باران‌سنجی خوب در نقاط بارش 
و استفاده از اطلاعات مکانی رادار برای بهبود ارزیابی مکانی بارش امکان‌پذیر است (اسلازینکا و زاویسلاک ۲۰۰۵ 
۳ باردوسی و داس ۲۰۰۸۵: ۸۵ بندینت " و همکاران۲۰۱۳: فصل اول ). 

۳1۳-5 مدل بارش "روانابی است که پارامترهای یکپارچه و همچنین پارامترهای توزیعی را پشتیبانی می‌کند 
(آهن ۲۰۰۷۵ ۱۷۹). در مدل یکیارچه پارامترهای هیدرولوژی به صورت متوسط بیان می‌شوند و تغییرات مکانی 
آن‌ها مد نظر نیست . یک مدل توزیعی به زیر حوضه‌هایی تقسیم می‌شود و پارامترهای هیدرولوژیکی برای هر واحد 
به طور مجزا تحلیل می‌شود. 

در تحقیقی که | سدی و بو ستانی (۲۰۱۳: ۱۱۹) به‌منظور برر سی عملکرد مدل ۳1۳-۳1315 انجام دادند. از مدل 
یکپارچه و مدل توزیعی استفاده کردند. همان‌طور که به صورت عمومی انتظار می‌رفت نتایج نشان داد که مقدار دبی 


اوج و حجم رواناب در مدل توزیعی عملکرد بهتری نسبت به مدل یکپارچه نشان داده است. 
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یکی از قابلیت‌های مدل‌هایی نظیر ۳1۳0-1۷15 برر سی ارات تغییرها ست. برخی تغییرات مکانی و برخی زمانی 
هستند. امروزه به دلیل افزایش جمعیت. تغیر کاربری ارا ضی با پو شش طبیعی به مناطق شهری به صورتقابل توجهی 
افزایش یافته بنابراین ایجاد این تغییرات و به دنبال آن کاهش پو شش گیاهی و تبخیر و تعرق در اندازه رواناب حاصل 
(اولانگ وفورست ۲۰۱۱ ۸ . دراین‌ارتباطمراد نژادی و همکاران (۱۳۹۶: ۳۵( در رابطه با الویت بندی زیر 
حوضه‌های مناسب بهره‌برداری جنگل با استفاده از مدل 1۳00-11۷15 آبه مطالعه پرداختند. مکان‌یابی مناطق مناسب 
جهت اجرای عملیات بهره‌برداری نشان داد به ترتیب زير در چه الویتی فراردارند. 
کاهش دبی اوج رواناب و جاری شدن سیل در بالادست رودخانه سدار را در به کار بردن ترکیپی از پروژه‌های 
افزایش نفوذ یا ذخیره رواناب اضافی دانست که با استفاده از مدل توزیعی ۲1۳-۳65 این مطالعه صورت پذیرفت. 
(الکساندر 7 ۲۰۰۹: ۱۱۸). 

سلیمانی و همکاران (۱۳۸۷: 6۸ درمطالعه‌ای با استفاده ازمدل ۱۳-5 آبه شناسایی مناطق سیل خیز پرداختند 
و به این نتیجه رسیدند که میزان مشارکت حوزه‌ها در سیل خروجی متناسب با دبی اوج نبوده ولزوما زير حوزه‌هایی 
پارامترهای مدل به‌طورعمومی * بهینه شده با شند. به‌عبارت‌دیگر به‌منظور دریافت پاسخ مناسب از مدل هیدرولوژیکی 
پارامترها باید ب‌صورت چند معیاره کالیبره شوند. مدل‌ها برای مناطقی با شرایط متفاوت با شرایط آب‌وهوایی و 
قنایو گرآفی مقتارک ازافه شله‌اند به همین حل باید انا کی شوه اما اکر کی اون بهضورت وی ورن 
نگیرد چه اتفاقی می‌افتد. هدف اول این تحقیق پاسخ علمی به این سوال است. 

اکثر تحقیقات صورت گرفته در زمینه کالیبراسیون مدل‌های هیدرولوژیکی با پارامترهای محلی است ولی اگر 


پارامتر ها در مقیاس های بز رگ‌تر کالیبره شو ند (تروی" و هم‌کاران۲۰۰۸: ۸ پوخرل و گویتا۲۰۱۰: ۱۲ 
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بلوسچل ۲۰۰۱۰: ٩۷۱۰‏ عبدالّه و لتنمایر ۱۹۹۷: ۲۵۳) می‌توان نتایج را به مناطق فاقد آمار تعمیم داد. برحی محققان در 
این زمینه تحقیقاتی انجام داده‌اند از قبیل سامانیگو (۲۰۱۰: ۲۰) پوخرل (۲۰۰۸: ۱۱ ) به این طریق که پس از مشخص 
کردن پارامترهای حساس؛ پارامترها به‌صورت منطقه‌ای کالیبره می‌شوند (هانداچا و باردوسی ۲۰۰۶: ۲۹۳؛ پاراجکا و 
همکاران۲۰۰۵: ۱1۸؛ بون ۲۰۰۱: )٩‏ و براساس پارامترهای محلیء پارامترهای بزرگقیاس را به دست آورده و 
کالیبره شوند. هدف از این مطالعه مقایسه واسنجی محلی و بزرگ مقیاس است. با توجه به اینکه در ایران کمتر به این 


موضوع اهمیت داده شده می‌توان نقطه قوت این تحقیق دانست و نتایج اين مطالعه را به مناطق فاقد آمار تعمیم داد. 


ای روف ها 

مشخصات حوضه موردمطالعه 

این مطالعه در حوضه آبخیز کوه سوخته که بخشی از حوضه آبخیز بهشت‌آباد است. انجام شده است. بخش‌هایی 
از شهرستان‌های شهر کرد بروجن و کیار در این محدوده قرار می‌گیرد. از نظر جغرافیایی بین ۰ درحه و ۶۰ دقبقه تا 
۱ درجه و ۲۰ دقیقه طول جغرافیایی شرقی تا ۳۱ درجه و ۲۰ دقیقه و ۳۲ درجه عرض جغرافیایی شمالی واقع شده 


مترومتوسط شیب منطقه /۱٩‏ است. همچنین میانگین بارش سالیانه در ایستگاه سینویتیک شهرکرد ۳۲۰ میلی‌متر است. 
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روش تحقیق 

در مرحله اولیه آماده سازی داده‌ها با استفاده از الحاقیه 1۳0-01301115 آبر اساس مدل رقومی ارتفاع ((۳) و 
گذراندن مراحل متفاوت. مشسخصات فیزیکی حوضه تهیه و به مدل 1150-1۷15 معرفی شد. در این مطالعه از 
حصوصیات فیزیکی و دینامیکی حوضه استفاده گردید. پارامترهای مدل به‌ صورت محلی و همچنین با در نظرگرفتن 
مقدار ثابت برای هر زير حوضه در حالت بزرگمقیاس. کالیبره شسدند. پارامترهای مربوط به روش تلفات» روش 
انتقال و روند یاپی جریان محاسبه و وارد مدل گردید. در مدل هواشناسی و مدل کنترل سازوکارهای لازم برای اجرای 
صورت گرفت در بخش سری‌های زمانی این مدل داده‌های بارش + ساعته و دبی * ساعته استفاده‌شده است. انتخاب 
داده‌های 1 ساعته به این دلیل بود که مطابق با رابطه تجربی اشنایدر (مدت زمان بارش) زمان بارش باید مساوی یا 
بیش از زمان هیدروگراف واحد با شد و به دلیل در د سترس نبودن داده‌های یک ساعته بارش, از داده‌های " ساعته 
| ستفاده گردید. به‌منظور شناخت پارامترهای < ساس تحلیل سا سیت پارامترها انجام شد. سپس برای برر سی 
و ضعیت کارایی مدل با حالات مختلف.1 رویداد برای واسنجی و دو رویداد برای صحت سنجی مدل انتخاب شد. 
روند مراحل محاسبه پارامترهای مدل 1۳0-11015به شرح ذیل است: 

مدل حوضه 

در این قسمت به‌منظور محاسبه پارامترهای موردنیاز در مدل تلفات و مدل انتقال در زیر حوضه‌ها به ترتیب از 
روش شماره منحنی 505 و هیدروگراف واحد ۹05 استفاده شد. پارامترهای موردنیاز مدل تلفات شامل نگهدا شت 
اولیه شماره منحنی و نفوذناپذیری است و پارامتر موردنیاز مدل انتقال زمان تأخیر است. 


با استفاده از روش برانسبای- ویلیامز زمان تأخیر محاسبه گردید (علیزاده ۱۳۹۰: ۱۹۰). 


۱ 09612 0) 
۲ 2 02 1 


۱۱۳۰ 1 (۲( 


که در این معادله. 16 زمان تمرکز برحسب ساعت. 1 طول آبراهه اصلی برحسب کیلومتر 13 اختلاف ارتفاع 
برحسب مت ۸ مساحت حوزه برحسب کیلومترمربع» یز" زمان تأخیر برحسب ساعت است. 

به صورت تئوری مقدار آ) بين صفر تا ۱۰۰ متغیر است اما در عمل بین ۳۰ تا ۱۰۰ می‌با شد. درآ برابر صفر 
روانابی از بارندگی حاصل نیامده و در آبرابر 1۰۰ تمامی بارش در سطح زمین جریان یافته و ارتفاع رواناب برابر 
ارتفاع بارندگی خواهد بود (مهدوی,۱۳۸۱: ۱۶۰). با استفاده از نقشه شماره منحنی که نشان‌دهنده خصوصیات حوضه 


ازنظر نفوذپذیری است. عدد مربوط به هر زیر حوضه با استفاده از ميانگین وزنی محاسبه شد. 


33۳5 بررسی اثر خصوصیات مکانی حوزه آبخیز بر واسنجی وقایع ... 


٩‏ حداکثر توان نگهداشت مربوط به نفوذ در خاک و ذخیره سطحی است. مقدار ٩‏ در رابطه با نوع پوشش و 
نحوه بهره‌برداری از اراضی و وضعیت سطح خاک ازنظر نفوذپذیری و داخل خاک ازنظر انتقال است. چنانچه ارتفاع 


نظر گرفته می‌شود (مهدوی»۱۳۸: ۱۶۸). 
۳ نم ۱۳۳ 
3 401 - 25400 
0 


در رابطه فوق مآنگهداشت اولیه و نگهداشت کل آال) شماره منحنی است. 
در مدل 100-۳1۷5 )11 آبعد از ایجاد مدل هوا شنا سی باید روش ارائه بارش م شخص شود که اینجا به دلیل 


موجود بودن هایتوگراف در این مدل هیتو گراف مشخص شده توسط کاربر انتخاب گردید. 
منظور مقدار این دو پارامتر از ۰- تا ۶۰ با فواصل ۰ تغییر داده شد و اثر آن بر دبی اوج سیل اعمال گردید. 


۱۱6۲۷۷7 - 0 (6) 
٩۳۱ - 10 
010: ۳6 


در رابطه فوق ٩۳۳۲‏ مقدار حساسیت تابع هدف حاصل از تغییر در اندازه پارامتر ۱۵ مقدار دید خروجی 
مدل حاصل از پارامتر جدید. 00 خروجی اولیه مدل. ۳ قدر مطلق درصد تغیرت است. 

واسنجی 
با توجه به اینکه دقت پارامترهای ورودی شاخحصی برای میزان کارایی مدل‌سازی است و نتایج شبیه‌سازی نهایی 
مدل به این مرحله وابسته است. بعد از تعیین حساسیت مدل ذسبت به پارامترها؛ پارامتر حساس انتخاب گردید و 


پارامترها به صورت محلی و بزرگ‌مقیاس کالیبره شد. در این پژوهش از پارامترهای نگهدا شت اولیه نفوذناپذیری و 
زمان تأخیر و همچنین"رویداد بارش برای واسنجی مدل استفاده شد. 

صحت سنحی 

ازانجاکه نایج حا صل از مدل در صمیمگیری‌ها و طراحی منابع آب‌وخاک و همچنین مسائل مربوط به سیل و 
امثال آن کاربرد فراوانی دار همواره درجه اعتبار آن‌ها مورد سوال است. با عدم قطعیت‌ها و خطاهای مدل‌ها از هیچ 


( ۷ 


۷ جغرافیا و مخاطرات محیطی شماره بیست و سوم 


مدل کامپیوتری نمی‌توان پیش‌بینی کامل و دقیقی را انتظار داشت. لذا آن‌ها را باید کالیبره و اعتبار سنجی نمود. صحت 
سنجی با یک سری اطلاعات جدید. مستقل از اطلاعات استفاده‌شده در مرحله واسنجی شبیه‌سازی جریان انجام شد 
که برای این پروژه صحت سنجی مدل از طریق اجرای مدل با پارامترهای بهینه شده برای دو رویداد بارش صورت 
گرفت. 

ارزیایی کارایی مدل 

برای ارزیایی عملکرد از شاحص ۹۷10111015 -«اععآبه‌عنوان معیاری برای سنجش تفاوت میان داده‌های مشاهده‌ای 
و شبیه سازی شده است هرچه مقدار شاخص به یک نزدیک‌تر باشد داده‌های شبیه‌سازی شده به داده‌های مشاهده‌ای 
شباهت بیشتری دارند. 


0( 0-۳ 
۴ 10-0 
در روابط فوق تا شاخحص ع1111ماناصمعلل 1 داده‌های مشاهده‌ای, ۳ داده‌های شبیه سازی شده ۵ا)متو سط 


5 -< 1 -- 


داده‌های مشاهده‌ای است. 


۳ بحث و نتایج 
نقشه مدل حوضه با استفاده از الحاقیه 11۳0-060-۳115 حاصل گردید (تصویر ۲) و نقشه حاصله به‌عنوان اولین 


0 260000 270000 280000 ۰ ۰-2909000۰ ۰-300000 ۰-3310000۰ 320000 
1 0 2 1 1 9 1 لس 


/ / 


/ 


3 

1880000 7۳ 890000 7 900000 7*۳ 910000 7*1 920000 7*7 930000 77 0000 77 0000 
1 
1 

1880000 ۰ 890000 900000 "910000 7" 920000 "930000 ۰ 000 7 00 


۲ ۲ 7 سسم و بو سس یر موی 
250000 260000 270000 20000 290000 23200000 310000 220000 


تصویر ۲. نقشه شماتیک حوضه کوه سوخته به‌منظور ورود به مدل ۲1۳0-11۷15 


بررسی اثر خصوصیات مکانی حوزه آبخیز بر واسنجی وقایع ... 
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20000 
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راهنما 


خطوط آبر اهه 
غلولانی ترین آبراهه 
زیرحوزه ۲ .. . | 
زیرحوزه ۴ ۱255 
زیرحوزه۱ [ 7 
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1880000 7 890000" 900000777 910000 7*۳ 920000 ۳ 930000 77 940000 ۳ 00 
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-320000 ۰-3310000۰ ۰ 300000 ۰ ۰.0 


تصوير ۳. زیرحوضه‌های حوضه کوه سوخته 


1 ۰۳ ۶2 ۳ 


۷۵ 


تحلیل حساسیت بر روی پارامترهای مدل نشان می‌دهد که پارامتر نگهداشت اولیه و زمان تأخیر دارای حساسیت 


کوچک در پارامتر موردنظر باعث تغییرات زیادی در مدل می شود و بالعکس اگر شیب نمودار کم باشد تغبیرات زیاد 


)۵ دشان داده شده مدل در هر دو حالت محلی و بزرگ‌مقیاس دسبت به پارامتر زمان تأخیر حساسیت بیشتری برای 


همه وقایع نشان می‌دهد. 


0 
100- 50- 0 
درصد تعسیرات 
]سوب اه 


محلی حوضه کوه سو خته رویداد ۰«۸(«۰(«۰(۳/۸« 


50 


1۳ ۳ ۳ ۳ 
100 100- 50 لا 50 
وف ت 

درصد تعسیرات 


عج|آسو واسهب 


دیی اوچ 


100 
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تصویر ۵. تحلیل حساسیت پارامترها در کالیبراسیون 


بزرگ‌مقیاس حوضه کوه سوخته. رویداد ٩۰/۰۵/۲۸‏ 


۷۳۹ 


جغرافیا و مخاطرات محیطی 


شماره بیست و سوم 


صاویر ذیل ذشان‌دهنده دبی اوج مشاهده‌ای و دبی اوج شبیه سازی شله است. صاویر 14 و ۰( مربوط به 


مقایسه دبی اوج مشاهده‌ای و شبه‌سازی شده در کالیبراسیون با استفاده از پارامترهای محلی تست و همچنین تصاویر 


(۷۵و 05 اختلاف بین دبی ارج مشاهده‌ای و شبه‌سازی شده در کالیبراسیون بااستفاده از پارامترهای بزرگ‌مقیاس را 


0 12 
۱ نا ٌ 
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۳ سس ۱ 
5 0 8 
3 : 
1057 
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تصویر ۷ دیی مشاهداتی و شبیه‌سازی شده در 
رویداد۰/۰۸/۲۹٩‏ در حوضه کوه سوخته-زیر حوضه 


(( کالیبراسیون بزرگ‌مقیاس) 


0 
قاگنا ۳ 4 ۳9 
۱ 1( 
51۳ 4 
د ۳ 
2 میت ۵ وکا | 3 
ِ ۱ 1 
بت 1 ث‌‌ِ 
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۱ 0 دا 
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8 .۰ 12:00 712 224 
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تصویر ٩‏ دبی مشاهداتی و شبیه‌سازی‌شده در 
رویداد ٩۱/۱۰/۰۱‏ در حوضه کوه سوخته-زیر حوضه 


(کالیبراسیون بزرگ‌مقیاس) 
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121۱/۵۱/۵۸۰ 
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تصویر آ دیی مشاهداتی و شبیه‌سازی شده در 


رویداد۰/۰۱۸/۲۹٩‏ در حوضه کوه سوخته-زیر حوضه 


۱( کالیبراسیون محلی) 


لا 


در 
‌ 
(1۵2101211)8000 
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تصویر ۸ دبی مشاهداتی و شبیه‌سازی‌شده در 
رویداد ٩۱/۱۰/۰۱‏ در حوضه کوه سوخته زیر حوضه 


۱( کالیبراسیون محلی) 


د 
لا 


0 


رویداد ٩۱/۱۰/۰۵‏ در حوضه کوه سوخته-زیر حوضه 


۱ کالیبراسیون بزرگ‌مقیاس) 


(حص)1211211 


لا 


دب 


19:12 


1124 


تصویر ۱ دبی مشاهداتی و شبیه‌سازی‌شده در 
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936 
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۱( کالیبراسیون محلی) 


1:48 


تصویر ۰ . دیی مشاهداتی و شبیه‌سازی‌شده در 


رویداد ٩۱/۱۰/۰۵‏ در حوضه کوه سوخته-زیر حوضه 


تصاویر فوق گویای برتری ظاهری کالیبراسیون محلی است. رویداد ۵ دی ٩۱‏ بیشترین دقت و رویداد ۱۱ 
رورش ٩۱‏ سیف وا کاسزه ما دک قیال این دیف مک اسر کی تابر ارام 
در فصل زمستان و بهار باشد. در کالیبراسیون بزرگ‌مقیاس نیز در رویداد ۵ دی ٩۱‏ کارایی بالایی مشاهده می‌شود. در 
کم انتون ند انیب از ات ارات هام مره بطان هه با راشای فان اف گیدافت ش 


اولیه نفوذناپذیری و شماره منحنی واسنجی گردیلن نتایج حاصل از واسنجی مطابق هل ( 2۷۱ یز امک, 


حدول ۱. پارامترهای بزرگ‌مقیاس واسنجی شده در حوضه آبخیز کوه سوخته 


زمان تأخیر (00) 
۰۸/۵ ۱۱۰ 
۱۰/۷ 


۲ 


تقوختانتیرس زان 
۱/۳۵ 
۲ 


/ 


نگهداشت اولیه (0۲0) 
8۰ 
۵۰ 


6۰ 


شماره منحنی ()) 
1۹ 
۸ 


1۹ 


مساحت (0000 
۳۵2 
۱۶۵۱/۵ 


۱۳۳ 


زیر حوضه 


پارامترهای موردنیاز اجرای مدل مطابق با آنچه در قسمت مواد و روش‌ها ذکر شد محاسبه گردید. به‌منظور دقت 


دک تن اولیه نفوذناپذیری و شماره منحنی مطابق اطلاعات جدول (۲) در "رویداد حوزه به‌صورت محلی 


واسنجی گردیدنك: 


۷/۸ جغرافیا و مخاطرات محیطی شماره بیست و سوم 


حدول ۲. پارامترهای محلی واسنجی شده11]:0-1115 در حوضه آبخیز کوه سوخته 


زمان تأخیر (0010) | نفوذناپذیری (/) ای ۳ مساحت (100) زیر حوضه رویداد 
۱ ((0) رصم 
۱۳۰ ۱/۸ 1۹ 5۵ ۱۳ ۱ 
۱۱6۰ /۳ 5۸ 0۰ ۱:۱۵ ۲ ۸/۹ ,۸« 
1۳ ۰/۹ نا 1۵ ۱۰۳۳/۹۲ ۳ 
۱۱۳۰ ۱/۵ 2۹ و ۹/۳ ۱ ۱ 
۱۰6۰ ۲ 0۸/۸ 0۰ ۱۱/۵ ۲ ۰/۱۱۱/۸۱۲« 
۹۸۰ ۱/۵ 11 0۲ ۱۰۳۳/۹۲ ِ 
۱۳۰ ۱/۲ 1۹ ۳۰ ۱/۳ ۱ 
۱۳6۰ ۱/۲ 11 ۳۰ ۱۵ ۲ ۱-۰۱۸۱ 
۱2۰۰ ۱/۵ 5۸ ۳۰ ۱۳ ۱ 
۱۰۰۰ ۰/۹ 24 5.۰ ۱۹/۳ ۱ 
۹.۰ ۱ 5۸ 0۰ ۱:۱۵ ۲ ۹۱۳۸*۳۹۸۱ 
۱۰۰۰ ۱ 11 0۰ ۱۰۳۳/۹ ۳ 
۳۰ ۱/۳ 1۹ ۳۰ ۱۳ ۱ 
۹۰ ۱/۸ ۸ ۳۰ ۱2۱/۵۵ ۲ ۱۳ 
1۹۰ ۱/۳ 1 ۳۰ ۱۰۳۳/۹۲ 3 
۱۳۰ ۱/۸ 1۹ ۳۰ ۱۳ ۱ 
۱۱۵۰ ۳/۸ 0۸۳/۸ ۳۰ ۱2۱/۵ ۲ ۹۱۹۱/۸۱/۵ 
۱۳۰۰ ۱/۹ 1 ۳۰ ۱۰۳۳/۹۲ ۳ 


نتایجم حاصل از واسنجی در * رویداد و حجم رواناب مشاهده‌ای و شبیه‌سازی‌شده در حوضه آبخیز کوه سوخته 
مطابق با جدول (۳) مورد مقایسه قرارگرفته است. 


حدول ۳ مقایسه نتایج واسنجی حجم رواناب به دو روش بزرگ‌مقیاس و محلی در مقابل مقادیر مشاهده‌ای 


مشاهده‌ای شبیه‌سازی 
حجم (۳) حجم () زیر حوضه رویداد 
واسنجی محلی واسنجی بزرگ‌مقیاس 
۹0 ۰/۸۱ ۰/۱ ۱ 
۰,۸۹۹« 
0 ۰/۷۹ ۰/۷۳ ۲ 


سال شه بررسی اثر خصوصیات مکانی حوزه آبخیز بر واسنجی وقایع 5 


ادامه حدول ۳ 
مشاهده‌ای شبیه‌سازی 
3 حجم ( 
حجم (2۳) 
واسنجی محلی 
۰/۵ ۰/۸ 
۱,۱ ۱/۳۵ 
۱/۸ ۱/۳۳۸ 
۱/۸ ۱/۳ 
۰/۳ ۰/۹۵ 
۱/۳ ۰/۹ 
۱/۳ ۰/۳۹ 
/ ۳۸ 
/ ۰/۳۱ 
/ :/ 
/ ۳۳ 
/ ۰/۳ 
/ ۳ 
۰/۵ ۰/۵ 
۰/۵ ۰/۹ 
۰/۳۵ ۳۳ 


واسنجی بزرگ‌مقیاس 


۰/۷ 
۱/۳0۹ 
۱/۳۲ 


زیرحوضه 


رویداد 


۹*۱۳," 


(۱۳*۱۸۱ 


۹۱۳/۰۹۹/۱ 


۱۰۹/۱ 


(۱/۱/۳۵ 


۷۹ 


پس از واسنجی مدل در 1 رویداد در حوزه صحت سنجی انجام گرفت. دو رویداد به‌منظور صحت سنجی مدل 


انتخاب گردید. دبی رواناب و همچنین حجم رواناب برآورد گردید. در جدول (۳) مقادیر سیر 


و داده‌های مشاهده‌ای ثبت‌شده در ایستگاه‌های هیدرومتری جهت مقایسه ارائ‌شده است. 


جدول ۶. ارزیابی کارایی مدل با استفاده از کالیبراسیون محلی و بزرگ‌مقیاس حوضه کوه سوخته 


تاریخ 
۱ 


۲( 
او( 
۱( 
۵( 


(۵ 


مرحله 


کالیبراسیون محلی 
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ادامه جدول ۶ 
5۳( تاریخ مرحله 
۳۱ 1۱-۳-۰۹ 
۰/۸۶ ۱-۲( 
۳ ۱۳۱ کالییراسیون 
۳۳ ۱۳۱ بزرگ‌مقباس 
۰۳ ۱۱۳۸ 
۰/۸۶ ۱۳۵ 


نتایجم حاصل از واسنجی و بهینه‌سازی نشان داد که برآورد دبی اوج و زمان رسیدن به دبی اوج هیدروگراف» 
مدل‌سازی با استفاده کالیبراسیون محلی برتری ظاهری دارد (جدول ۳ و ۶)(تصاویر ۸۰۱ و ۱۰) هرچند این نتیجه 
برای کل حوضه عمومیت ندارد. همان‌طور که انتظار می‌رود از لحاظ ارتفاع رواناب کمترین اختلاف بین رواناب 
مشاهده شده و شبیه سازی شده مربوط به کالیبرا سیون با استفاده از پارامترهای محلی است (جدول ۲). در ارزیابی 
کارایی نتایج. شاحص 0ذ01ا50-(9هآلبرای کالیبرا سیون محلی ۰/۸۵ و کالیبرا سیون بزرگقیاس ۰/۱۵ به‌د ست‌آمده 


که در دامنه مطلوبی برای شبه‌سازی قرار دارد. 


۶ نتیجه گیری 

مدل سازی از روش‌های غیرهم‌ستقيم است که سرعت بیشتری نسبت به روش‌های میدانی دارد. برای دستیابی به 
نتایج دقیق در مدل‌سازی نیاز به برآورد قابل‌قبول پارامترهای مدل و نیز اطلاع از متغیرهای زمانی و مکانی اسست. به 
این فررآیند کالیبراسیون اطلاق می‌شود. نتایج حاصل از مدل سازی بارش-رواناب توسط مدل ۳1۳0-۳115 در حوزه 
موردمطالعه. گویای این مطلب است که استفاده کالیبرا سیون محلی نسبت به کالیبرا سیون بزرگقیاس دقت ظاهری 
بیشتری دارد و این موضوع مطابق با تحقیقات سامانیگو و همکاران (۲۰۱۰: ۲۲) و هاندچا و باردوسی (۲۰۰۶: ۲۹۳) 
است؛ اما مقدار شاخحص محاسباتی ۷26-00011616 در کالیبراسیون بزرگ‌مقیاس نسبت به کالیبراسیون محلی 
تغییرپذیری کمتری دارد. دبی اوج رواناب شبیه سازی شده در کالیبرا سیون محلی در مقایسه با رواناب مشاهده‌ای و 
همچنین حجم رواناب شبیه سازی شده در کالیبراسیون محلی در مقایسه با رواناب مشاهده‌ای اختلاف کمتری نسبت 
به این پارامترها در کالیبراسیون بزرگ‌مقیاس دارد که اين مقایسه نشان از برتری ظاهری کالیبراسیون محلی دارد. 
بااین‌وجود هر دو کالیبراسیون در برآورد رواناب پاسخ‌های قابل قبولی ارائه دادند در حالی که پارامترهای به‌دست‌آمده 
در حالت عمومی يا بزرگقیاس حساسیت مکانی کمتری دارند. این موضوع امکان تعمیم نتایج واسنجی را به مناطق 
فاقد آمار بیشستر می‌کند. لذا همان‌طور که در بخش دوم اين مطالعه از روش کالیبراسیون با استفاده از پارامترهای 
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بزرگ‌مقیاس نشان داده شد برای پروژه‌های میدانی توزیعی و پژوهشی مناسب می‌باشد ولی به‌طور موضعی کارایی 
روش محلی بیشتر بوده است که این نتایج پاراجکا (۲۰۰۵: ۱1۸) هماهنگ است. هر دو روش در حوضه کوه 
سوخته با مساحت یکسان استفاده گردید. نتایج نشان می‌دهد که این دو روش واسنجی با هم تفاوت قابل‌ملاحظه‌ای 
دارند و هزینه‌ای که پارامتره کردن بزرگقیاس استفاده می‌شود درواقع صرف کاهش حساسیت به وابستگی مکانی و 
ایستایی بهتر پارامترها می‌شود که منجر به تولید نتایج قابل‌اعتمادترمی‌شسود. به‌طورکلی با استفاده از پارامترهای 
بزرگ‌مقیاس در زمان صرفه‌جویی می‌شود که با نتایج تروی (۲۰۰۸: ۸ 4 پوخرل (۲۰۰۸:۱۱ ) و بون (۲۰۰۱: )٩‏ 
هم‌خوانی دارد. تحلیل سا سیت عوامل ذشان می‌دهد پارامترهایی مانند نفوذناپذیری و نگهدا شت اولیه حساسیت 
ب سیار بالاپی را در منطقه دارند. همچنین مدل ذسبت به تعیین زمان تأخی عمل واببستگی زمانی را ذشان می‌دهد که 
رسیدن به حالتی عمومی را دشسوار می‌کند. تحقیقات تکمیلی برای رسیدن به زمان تأخیر قابل‌اعتماد به صورت 


بزرگ‌مقیاس برای شبیه‌سازی رواناب برای مطالعات بعدی پيشنهاد می‌شود. 
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9)3(,137-۰ ,56060195 5661665 551۵1۱ 60۳1۲ ۵ و۲۵۱0 6۲5۰)عصصه تهج معط عحمحصطم‌اهع 

۰6 ها مع1عماه0او م220ه نع لهتاومو ۵۶ عون معط 0۵ ,(2010) ۷۰ ۳۱ مقانتات) عک ر.ظ بلعتط0ظ 
46)1(۰ ,6۹60۲61 ۳۵۵۵۵۰۵۵5 ۲۷۵۲۵ ۱0۵0615۰ ۷۷۵۵۵۲۹۳۴۵0 015010۵0 ۵۶ ط0تاه1ا1هع 

منهج ما طمه200(۲0 ۲6۵۱2۲1220100 [0212: ۸ 2006(۰) .1۰ ,۷۷۵۵6۵0۵61 ی ,۷۰ .۳۱ ما۷0 ر.ظ بلعت0ظ 
4)12(۰ 6۵۵۵۲6۱4 ۵90۱۵۵۵5 ۷۷۵/۵ 80061۰ ۳/2۵90 60)اهاتتاوتل 2 ۲۵۲ ممتاقصصتاوع 

-110ع ۵ 01 ۱6۵100۵112210100 تمامصه نه0 ۱۷۲۵۱0۵6۵18 2010(۰) 5۰ رتععط1اا » ریک رتناک رم ,52102016820 
1-۰ ,(46)5 63601 ۲۵509۲۵۵5 ۷۷۵۶۵۲ ۱0۵9056۵16۰ 66 21 ۱۱۵061 0۳010216 2960 

مافصتاه 10۲ فعمتامحعهاصا متمطامومصاه-صه1 ز۵ ما۵ م1 .(2008) ب ماه عک ویو ممص0860۷1۲۵0 
ب(33)1 ,۲6۵۲5 100 عوظ ۱۲6 اوعد که ۷۵۱۵۵ ۱۸۵۸6 ,۲۱۲00۵ ها رازه۱۷2 
179-4 

۲ ۲۱۱۷۱۵-)۲۱۳ ومتودا عصتام0 مج ۶۴مصیاااحتصتقط ,(2015) .۱۲ ۷۲۰ ,عافک عک رب۷۲ .9 ,عله] ریک ۸۲ ,نا 
,1661110108165 01۱8 616166 ۱۲18,۹ 61811۱66۳ 1 ۲۵۹۵۵۲6۲ 0۴ ۵۵۵۵ ۱۵۳۱۵۵0۵ 02618۰ 1۱۷۵۲ 2 ۱۷۷ 
28)1(,21-0 

-لقاصمصتاجهع ۶۵۲ 0مطامجط مماهتمانلهه اصملترله مه ,(2008) .1 ملام1مط9 عک و.۲ .۲ ,۷۷۵۵0 و.؟ و1۲0۷[ 
.(44)9 ۵560۲61 ۵509۵۵5 ۲۷۷۵۵۲ ,3۱0001102 ومهگتناد 1200 86216 

۲ ۵7 50۱ .015002186 م0651 10۲ وام۵0و ۶مصااحگصنه ۵۶ م0ولنهم‌جصمن 2013(۰) .1 2890۷ ۲۷ 
26-۰ و(8)01 ,06560۲61 

1 ۲2۱0۶۵11 لمح ممناعناتاد م۵0 لهمنع1010 بط ۵ 0۶ 6۴۲6۵۲ 1۳6 .2015(۰) .با رکام16122 ک ر.ظ ,۷۷21682 
305-۰ ,(47)4 500۲۷ - 5667665 ]0 ۱۲۵ 0۱۵۵0۲ د۸۱۲۵ 

05 ۵۵21 م1600۵18 10۶ وملععاه0و طمتاهوتا له 1۵00 0۶ اممجطروعووه ,(2014) .ط) ,ع)لهءنا ما2 ۱۷۷ 

10۷2 ,10772 0۶ ۲ازوته۲ نصا (فتعمطا و عمافمصه 0مصوتامنام‌هنا) 0مصفماه ۳۳۷۵ بقلم موروزنا فط رز 

( 


۸۶ 


جغرافیا و مخاطرات محیطی 


شماره بیست و سوم 


